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ABSTRACT 

We present 'BYTe', a comprehensive 'hot' line list for the ro-vibrational transitions 
of ammonia, 14 NH 3 , in its ground electronic state. This line list has been computed 
variationally using the program suite TROVE, a new spectroscopically-determined 
potential energy surface and an ab initio dipole moment surface. BYTe, is designed 
to be used at all temperatures up to 1 500 K. It comprises 1 137 650 964 transitions in 
the frequency range from to 12 000 cm^ 1 , constructed from 1 366 519 energy levels 
below 18 000 cm^ 1 having J values ^ 36. Comparisons with laboratory data confirm 
the accuracy of the line list which is suitable for modelling a variety of astrophysical 
problems including the atmospheres of extrasolar planets and brown dwarfs. 
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1 INTRODUCTION 



Ammonia is the main nitrogen-containing molecule observable in a variety of astrophysical environments. For example, whilst 
the spectra of M and L-type brown dwarfs are dominated by H2O, and CH4 becomes increasingly important at lower tem- 
peratures, NH3 is a significant source of opacity in the atmospheres of late T-type dwarfs, particularly in the 10.5 /im region 



(Burrows et al. 1997 Sharp & Burrows|2007). Modelling suggests that absorption by ammonia will be even more important 



in the yet-to-be discovered Y-dwarfs ( jBurrows et al.|2003 1 and it is possible that NH3 absorption bands may be the principal 
distinguishing feature of this new class of ultra-cool dwarf. 

The molecule is known to be present in the atmosphere of solar system gas giant planets ( Lara et al.|1998 |. For example, 
emission spectra of hot ammonia were observed in the atmosphere of Jupiter after the impact of Comet Shoemaker-Levy 9 
( Orton et al.|1995 I and in the aftermath of the July 2009 impact on Jupiter ( Orton et al.|2010 Fletcher et al.||2010 1 

NH3 is present in cometary comae, typically in number densities ~ 0.5% of that of gaseous H2O (Kawakitafe Watanabe 
J2OO2 



Mumma 



Boncv et al. 2009). Although the cometary environment is cold, solar pumping IjDello Russo et al.|2004 



Kawakita & 



2010) and other excitation mechanisms ( Barber et al.||2009 l are able to excite high lying ro-vibrational molecular 



states that in conditions of thermal equilibrium would only be significantly populated at temperatures in the region of 2 000 
K or higher. In the densities that exist in the coma, these states are able to decay radiatively before being collisionally de- 
excited, and spectral lines originating from these highly excited states give important clues to the physical processes at work 
in cometary comae. The Einstein A coefficient listed in BYTe will enable the measured intensities of cometary NH3 spectral 
lines to be used in modelling the conditions in the coma. 



Ammonia masers are also observed. These also involve transitions between high lying states (Madden et al. 1986) 



Modelling maser action requires significant spectroscopic data, and this too can be provided by our line list, which we call 
BYTe. 

An exciting new area of astronomy is the characterisation of extrasolar planets. Water and methane have already been 



shown to be present in exoplanet atmospheres (for example in HD 189733b (Tinetti et al. 2007 Swain et al. 2008 1), but 
to-date NH3 has not. The reactions by which N2 and H2 are converted into ammonia in the atmospheres of expoplanets and 
brown dwarfs are complex and outside the scope of this paper. However, the equilibrium between N2 and NH3 favours NH3 
at lower temperatures as (to a lesser extent) do higher pressures. These temperature and pressure dependencies suggest that 
the outer atmospheres of extrasolar giant planets at large orbital distances will contain significant quantities of ammonia 
( Sudarsky et al.|2003 l, and it is likely that NH3 will soon be confirmed to be present in these objects. 

Before BYTe, there was no NH3 line list that was sufficiently complete and accurate for use in modelling the atmospheres 
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of exoplanets, brown dwarfs or other astronomical objects at elevated temperatures. Many astronomers have used the experi- 
mental NH3 line list of ] 



Irwin et al. 



( 1999 1, which covers the 4 000-11 000 cm -1 spectral region, even though this line list is only 
designed to be used at temperatures below 300 K. The HITRAN database ( Rothman et al.||2009 l has been also extensively 
used. This contains approximately 30 000 lines of 14 NH3, but it too is only suitable for use at ambient temperature, or below, 
and even then is lacking data in significant regions of the spectrum; we recently demonstrated this using our computed, 
comprehensive line list for cold NH3 ( |Yurchenko et al.|2009| ). This line list comprises 3.25 million transitions between 184 400 
levels. It has an upper energy cut-off of 12 000 cm -1 and a maximum rotational quantum number J = 20. However, despite 
having many more lines than any other NH3 list, like those other lists, it was designed for use at temperatures up to 300 K, 
which renders it unsuitable for most astronomical applications. This problem is overcome by BYTe. 



2 THE BYTE LINE LIST 



BYTe is a catalogue of transitions represented by frequencies, Einstein coefficients, energy levels, and quantum numbers, 
which together fully characterise the electric dipole transitions of the ammonium molecule in the frequency range from 
to 12000 cm -1 . The list was computed variationally using the program suite TROVE ( Yurchenko et aL||2007 |. It comprises 
1 137650 964 transitions between 1366519 energy levels, all below 18000 cm -1 , and for all Js up to J = 36. All transitions 
are within the ground electronic state of the 14 NH3, since the excited electronic states are above the dissociation limit of the 
molecule. Less than 30 000 NH3 lines are known experimentally: BYTe contains about 40 000 times as many. It is the most 
comprehensive and accurate line list of NH3 and is designed to be applicable at temperatures up to 1 500 K. 

The ro- vibrational energies and wave functions were computed using a new 'spectroscopic' potential energy surface (PES) 
of NH3 (Yurchenko et al. 2011 1. This PES was recently generated through a least-squares fit to the 14 NH3 experimental energy 
levels mostly taken from the HITRAN database (|Rothman et al.||2009|). In order to evaluate the Einstein coefficients an ab 

is used. The reported 'hot' NH3 line list is the 



initio (ATVZ) dipole moment surface of NH3 from 



Yurchenko et al. 



(|2005J 



result of more than three years work by our team. The method of calculation was similar to, but more computer-intensive 
than in our earlier T = 300 K NH 3 line list ( | Yurchenko et al.|2009 |. 

Preliminary versions of BYTe have already been used in several astrophysical studies. In particular the data has been 
used in analysing the atmosphere of the late T dwarf UGPS 0722-05 ( |Lucas et al.|2010| ) and the atmosphere of the transiting 
hot Neptune GJ436b (|Beaulieu et al.|2010| 



3 COMPUTATIONAL DETAILS 



The variational rotation-vibration program suite TROVE ( | Yurchenko et al.| [2007) was employed for all nuclear motion cal- 
culations required to produce the BYTe line list. This procedure was the same as that used in our cold line list ( |Yurchenko| 
et al.|p009 l, and the reader is referred to this earlier paper for a more detailed account of TROVE. The main difference in 
generating BYTe was that a larger basis set was used than in the case of the cold NH3 paper. 

TROVE normally requires a polyad number P to control the basis set. In the case of ammonia, we define: 

P = 2(m +n 2 + n 3 ) + n 4 + n 5 + — , (1) 

where m are the quantum numbers associated with the basis functions <j> ni , see Yurchenko et al. (20091. Thus, we include 
in the basis set only those functions <j> n for which P ^ P m ax- Pmax was set at 14. This optimum value produced converged 
eigenvalues with the minimum computational demands. 



The Hamiltonian matrices were constructed from the ( J=0)-contracted basis set (Yurchenk o et al.|2009 |. In order to im- 
prove the agreement of our line frequencies with experiment, the empirical basis set correction (EBSC) was utilized ( |Yurchenko| 
et al. 2009). Using the EBSC approach, the vibrational energies in the ro- vibrational calculations were replaced by the corre- 
sponding experimental band centers, where these are available. 

Apart from extending the basis set, we also enlarged (i) the energy range from -E ma x = 12 000 cm -1 to iJmax = 18 000 cm -1 , 
(ii) the frequency range from - 8 000 cm -1 to - 12 000 cm -1 , and (iii) the range of the rotational excitations considered 
from J max = 20 to J max = 41, respectively. These changes were necessary for the new line list to be sufficiently complete and 
accurate to enable spectral calculations at temperatures up to T = 1500 K, which is far more demanding than at T = 300 K, 



which was the temperature criterion adopted in Yurchenko et al. ( 2009 ) 



The extensions to these key parameters has a huge impact on the computational costs. In the case of BYTe, the largest 
Hamiltonian matrix to be diagonalized was for J = 36 (E symmetry block) with dimension 162 763. In contrast, the largest 
Hamiltonian in Yurchenko et al. (20091 was 86 000 at J = 20. However, the most elaborate calculations, were in the J = 
20-27 range. These were characterised not only by large matrix dimensions (80 000 - 110 000), but also by a large number of 
roots of the corresponding eigenvalue problems. This is illustrated on Fig. [l] where we show the number of eigenvalues below 
the 18 000 cm" 1 threshold and the dimensions of the i?-symmetry matrices for J = 0. . .41. Using the computer facilities 
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available to us we could efficiently employ the LAPACK diagonalizer dsyev for all J below 23. For J = 24 and higher we had 
to switch to the Lanzcos-based iterative diagonalizer PARPACK ( |Lehoucq et al.fl998 l (a parallelized version of ARPACK). 
In fact, the dsyev LAPACK routine allowed us to compute and store all eigenvalues of a given Hamiltonian matrix for J ^ 23, 
beyond the -E max = 18 000 cm -1 threshold. One can see from Fig.[l]that for J > 30 the number roots to be determined in the 
eigenvalue problem represents only a small fraction of the matrix dimension (less than 10 %). This is important for efficient 
utilization of the iterative Lanzcos-type diagonalization approaches. The matrix diagonalization is the most consuming part 
of line list calculation, both in terms of time and memory. 

Using these techniques, we were able to generate all eigenvalues and eigenvectors for values of J up to 41, subject to the 
energy threshold of 18 000 cm -1 . However, only levels below J = 36 were taken into the line list because higher J values did 
not have any ro- vibrational states with energies below the threshold for the lower energy states in BYTe of E\ ovl = 8 000 cm - . 

As previously mentioned, in order to compute the eigenvalues and eigenvectors of NH3 we employed a new, improved, 
'spectroscopic' PES, NH3-2010. This PES was recently generated through fitting to the available experimental data below 
J = 8. It involved adjustments to the analytical representation from Yurchenko et al. (2005), using the refined PES from this 
work as a starting point. This refinement was performed using TROVE, which was extended to allow such fitting tasks. The 
details of the refinement procedure and the new PES NH3-2010 will be given elsewhere ( Yurchenko et al.|2011 l. 

An accurate dipole moment surface (DMS) is a prerequisite for producing accurate line intensities. Tests have shown that 
these are best taken directly from high quality ab initio calculations (Lynas-Gray et al. 19951. For BYTe we employed the 
ATZfc DMS of NH3 from I Yurchenko et al.| ( |2005[ ) in the improved DMS representation from | Yurchenko et al.| ( |2009| >, which 
behaves well when the molecule adopts a planar geometry. This DMS has been shown to provide intensities in good agreement 
with experiment (Yurchenko et al. 2005, 2009); it was also used in producing the 'cool' ammonia line list (Yurchenko et al. 
[2009] ). 

To speed up the computation of the transition moments used to generate the Einstein A coefficients, we used the pre- 
screening procedure, see Yurchenko et al. (20091, in which only eigen-coefficients with magnitude larger than 10~ 12 were 
selected. This reduced the size of the vectors (by about 70-80%) as well as the computation time. We also imposed an 
absorption intensity threshold 10 -12 cm/ mol (1.7xl0 -36 cm/molecule) at T = 1500 K, for the intensities to be included in 
BYTe. This threshold corresponds to about 10 -16 of the maximum intensity at T = 1500 K (18 000 cm/mol). The major 
bottle-neck in these calculations was associated with the need to read repeatedly a huge number of eigenvectors stored on the 
disk. The large size of these vectors prevented us from keeping all of them in the virtual memory during the computational 
process. We were able to reduce the number of readings by optimizing the computational logistics. In intensity simulations 
each transition can be independently processed, which effectively naturally parallelizes them, distributing them between 
computational nodes (see Yurchenko et al. ( 2009 1). 

The most expensive part of the intensity calculations was for J from 10 to 20, which accounts for approximately 75 % 
of all transitions to be computed. In Fig. [2] we show how the number of lines and energy levels depend on J. The number 
of levels peaks at J = 15 and the reduces gradually with higher J and eventually reaches zero due to our lower and upper 
energy limits of 8 000 and 18 000 cm -1 respectively. 



4 STRUCTURE OF THE LINE LIST OF NH 3 



The BYTe 14 NH 3 line list contains 1.138 billion lines, which are the allowed transitions above a certain minimum intensity, 



between 1.367 million ro-vibrational levels. The structure of BYTe is similar to that of the BT2 water line list (Barber et al 



2006). BYTe comprises two files. One, the Energy file, holds the energies and ro-vibrational quantum numbers of all NH3 
states up to J ^ 41 that are less than 18 000 cm -1 above the zero point energy (our value is 7430.288276 cm -1 ). Table[l]gives 
an extract from the Energy file. Apart from the general quantum numbers associated with the molecular group symmetry 
r ( Bunker fc Jensen|1998[ ) (NH3 belongs to X>3h(M)) and total angular momentum J, we have used both 'normal mode' and 
'local mode' quantum numbers when labelling the energy levels. Our labelling scheme is discussed below. 

The structure of the Transition file is simpler. It contains three columns: two give the reference numbers of the upper 
and lower states as they appear in the Energy file and third contains the Einstein A coefficient (s _1 ) for the transition (see 
Table |2J. The entries in the Transition files are sorted according to the frequency, and we have split them into 120 small files 
in order to reduce the amount of data that needs to be handled when examining a specific frequency region. 

In the actual TROVE calculations we employed 'local mode' basis functions in the FBR representation as explained in 



detail in Yurchenko et al. (20091. This allowed us to label the energy levels and hence to assign transitions based on the 



particular basis set making the largest contribution within the appropriate eigenfunction. Our local mode quantum numbers 
include F ro t , K, r rot , r v ib, ni, TI2, TI3, ri4, 715, ne- Here K is the projection of total angular momentum onto the molecular 
symmetry axis; m, 112, 713 are stretching local mode quantum numbers (Mills & Robiette 19851 which correlate with the 
normal mode notation as n\ + n% + 713 = V\ + 1/3; n,4 and 715 are deformational bending quanta; ng is the inversion quantum 
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Table 1. Extract from the BYTe Energy file. 
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Column Notation 



1 


N 


Level number (row) 


2 


J 


Rotational quantum number, angular momentum 


3 


r 


Total symmetry in (M) 


4 


^block 


Level number in a block 


5 


E 


Term value (in cm -1 ) 


6 ... 9 


nj . . . T14 


Normal mode vibrational quantum numbers 


10, 11 




Vibratiobal angular momenta, normal mode vibrational quantum 
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1"inv 


Inversional parity (0,1) 


13 


J 


Rotational quantum number (the same as column 2) 


14 
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Rotational quantum number, projection of J onto the z-axis 
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Trot 


R.otational parity (0,1) 


16 ... 21 


ui . . . v e 


Local mode vibrational quantum numbers (see (|Yurchenko et al.|2009^) 


22 


Tvib 


Symmetry of the vibrational contribution in 03^(1^!) 



Table 2. Extract from the BYTe Transition file. 
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number equivalent to 2v2 + n nv , where v 2 is the normal mode quantum number and n nv = ng mod 2 is the inversion 
parity ( Yurchenko et al.|2005 l. Finally, r rot and F v ib are the rotational and vibrational symmetries in Z>3h(M). 

Apart from the 'local mode' assignment that is generated by T ROVE, we also provide the standard normal mode quantum 
numbers Vi, V2, v 1 ^ , v 1 ^ , according with the Herzberg convention (Herzberg 19451. V\ and V2 are the symmetric stretch and 
symmetric bend, respectively, whilst V3 and «4 are the asymmetric stretch and asymmetric bend respectively. The additional 
quantum numbers £3 and h are necessary to resolve the degeneracy of the «3 and V4 vibrational states, respectively. 

The selection rules which determine the allowed electric dipole transitions of 14 NH3 are AJ = J' — J" = 0, ±1 (J" + J' ^ 1) 
with symmetry selection rules, A' 2 4-> A'2, and E' o E" . We used the nuclear spin statistical weight factor g ns = 12 and 6 
for the A' 2 <r> A'2 and E 1 o E" transitions, respectively. The A'\ and A" levels are characterized by g ns = 0, that is, the 
corresponding transitions do not exist. 

It should be noted that our assignments do not always agree with the experimental ones for the following reasons: (i) 
ambiguous definition of the quantum numbers (apart from J and T), which depend on the basis functions used; (ii) strong 
interactions between ro- vibrational states of close-lying levels; and (iii) mapping between the normal and local mode labels is 
not always straightforward. The last of these means that it is sometimes difficult to distinguish between the symmetric and 
asymmetric stretch quantum numbers V\ and ^3. 



5 LINE LIST: DESCRIPTION AND VALIDATION 

Absorption and emission intensity simulations are temperature-dependent. Specifically, temperature appears in the Boltzmann 
factors exp(—E/kT), where k is the Boltzmann constant. Line lists, in contrast, do not specify a temperature, since the 
Einstein coefficients for the transitions are independent of temperature. We therefore need to explain why we refer to BYTe 
as a 'hot' list. The reason is that BYTe is able to describe the absorption/emission processes in NH3 for temperatures up to 
Tm ax = 1 500 K. This is because the energy threshold that we have adopted, E\ OVJ — 8 000 cm -1 , ensures that we include all 
those states that are significantly populated up to T = T max . 

The temperature coverage of our data can be conveniently checked by using it to compute temperature-dependent partition 
functions for NH3. 
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Figure 1. Dimensions of the _E-symmetry matrices (squares), the corresponding number of eigenvalues below 18 000 cm 1 (circles) and 
number of non-zero elements on each row (diamonds). 



1x10 




Figure 2. Numbers of levels below 18 000 cm 1 per each J (lower display) and corresponding transitions for each J — > J + 1 (upper 
display). 



Q = Y,9i eM-Ej/kT), (2) 

i 

where gj is the total degeneracy of the state with energy Ej and the sum runs over all energy levels of the molecule. Including 
all energy levels up to 8 000 cm -1 gives an NH3 partition function, Q, of 57944 at T =1 500 K. This value is about 5 % below 
our best estimate for Q at 1 500 of 61 223. This discrepancy, which rapidly disappears at lower temperatures, is the reason we 
do not recommend the use of BYTe for temperatures above 1 500 K. 

Apart from E\ ovl there are two main contributions that affect the total number of the excited ro- vibrational states accessed 
in the variational calculations, (i) the size of the vibrational basis set, which BYTe is determined by the polyad number P max 
(see EqjTJ and (ii) the size of the rotational basis set, determined by the maximal total angular momentum J ma x- Fig. [3] 
illustrates how the partition function of NH 3 depends on J max for T = 300 K, 600 K, 900 K, 1200 K, and 1500 K (upper 
display). For T = 1 500 K the contribution to Q from J > 30 amounts to less than 1 %, see lower display of Fig. where the 
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Figure 3. The partition functions Qj max of NH3 at different temperatures vs. the maximum J value used in eq. |2jl, J ma x, (upper 
display) and corresponding contributions (Q — Qj max )/Q {%), where Q = Qj max — 41. 




Figure 4. Absorption spectra of NH3 given by BYTe for T = 300, 600, 900, 1200, and 1500 K, convoluted with a Gaussian profile, 
HWHM = 0.5 cm" 1 . 



relative contributions (Q — Qj max )/Q to the partition function are also shown as percentages. Here we used Q — Q./ max= 4i, 
where J ma x = 41, corresponding to the highest J in our variational calculations. 

The effect of the vibrational basis set on the partition function is less pronounced: at T = 1 500 K and J max = 35 the 
partition function is converged to better than 0.03 % when testing our choice of basis set size. There are, of course, other 
factors that give rise to errors in our computed partition functions. One of these is the potential energy function, which is 
only reasonably accurately fitted to low-lying ro-vibrational levels. We will address the nature of the partition function of 
ammonia in greater details elsewhere. 

Figs. [4] and [5] illustrate the T-dependence of the absorption spectra of ammonia, and give an overview of the complete 
range as well of the four selected regions computed at T = 300 K, 600 K, 900 K, 1 200 K, and 1 500 K. As expected, the 
spectrum profiles at higher T become less extreme as the populations from the lower vibrational states are reduced in favour 
of the vibrationally excited states. Only those few bands that are the most pronounced features at T — 300 K are still 
recognizable at T = 1 500 K. 
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Figure 5. BYTe absorption spectra of NH 3 for T = 300, 600, 900, 1200, and 1500 K, convoluted with a Gaussian profile, HWHM = 5 
em : The 3 fim, 2 /im, 1.5 (im and 1.25 fim regions. 
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Figure 6. Absorption spectra of NH 3 at T=300 K (cm/mol): BYTe vs HITRAN. 



Our 'cold' NH3 paper (Yurchenko et al. 2009 1 contained detailed comparisons with the HITRAN database (Rothman 
|et al.|[2009[ ) which demonstrated the accuracy of our procedure and also, that even at 300 K, HITRAN is missing significant 
NH3 data. In BYTe we have improved our ability to reproduce the HITRAN data, see Fig. [6] however we do not show the 
entire comparison here. 

As we were completing our calculations, a high temperature emission spectra of ammonia (T = 900 K) for the ground 



and V2 states of ammonia was reported by Yu et al. (2010 1. We use this spectrum to provide an independent validation of the 
BYTe line list. To this end we used BYTe to generate a 900 K synthetic emission spectrum of NH3 in the same spectral region. 
This is shown in Fig. [T] (lower part). The 'experimental' spectrum on this figure was generated using the Einstein coefficients 
from the synthetic line list reported by Yu et al. (20101. The spectra agree not only qualitatively but also quantitatively 



in terms of the absolute intensity values. The additional lines in the upper part of Fig. [7] are present because our synthetic 
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200 400 600 800 1000 1200 1400 1600 1800 2000 
wavenumber / cm -1 

Figure 7. Emission spectra of NH3 at T=900 K (Ergs/molecule/Sr). The 'experimental spectrum' (lower display) is reconstructed 
using the Einstein coefficients from |Yu et ah] | |2010| l . 



spectrum includes all possible transitions falling into the region including hot bands, while Yu et aFs data is only for the 
ground and 1/2 states of NH3. 



6 CONCLUSION 

We have calculated to a high level of accuracy the frequencies and Einstein A coefficients of all the transitions that are present 
in the emission and absorption spectrum of NH3. The only limitations are: upper states with energies above 18 000 cm -1 
are excluded; there is an effective short-end wavelength cut-off of 1 fj,m, due to the incompleteness of our data at frequencies 
above 10 000 cm -1 (the excluded region is unimportant for NH3), and extremely weak lines have been excluded, which is of 
little physical significance. 

Although the BYTe line list is explicitly aimed at modelling hot ammonia, it improves on our previous cold line list 
( Yurchenko et al.|2009 ) in terms of the quality of the potential energy surface used, in the size of the basis sets employed and 



in the range of frequencies studied. We therefore recommend the use of BYTe for all temperatures up to 1 500 K. The line list 
is freely available and can download in its entirety or in parts from: http://www.tampa.phys.ucl.ac.uk/ftp/astrodata/BYTe 
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